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ABSTRACT 



The structure and energy balance of the solar chromosphere remain poorly known. We have used the imaging spectrometer IBIS at 
the Dunn Solar Telescope to obtain fast-cadence, multi-wavelength profile sampling of Ha and Ca ii 854.2 nm over a sizable two- 
dimensional field of view encompassing quiet-Sun network. We provide a first inventory of how the quiet chromosphere appears in 
these two lines by comparing basic profile measurements in the form of image displays, temporal- average displays, time slices, and 
pixel-by-pixel correlations. We find that the two lines can be markedly dissimilar in their rendering of the chromosphere, but that, 
nevertheless, both show evidence of chromospheric heating, particularly in and around network: Ho- in its core width. Can 854.2 nm 
in its brightness. We discuss venues for improved modeling. 
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The solar chromosphere is one of the most difficult solar atmo- 
spheric regi mes to observe, model, and under s tand (see recent 
reviews by [Judge & Peteii 119981: lJudgd 120061: ICarlssonI l2007t 
lRuttenl200C In this paper we compare spectrally resolved chro- 
mospheric image sequences taken in the Balmer Ha line at 
656.3 nm and the Can 854.2 nm line (henceforth abbreviated 
to CaIR). Both lines sample the chromosphere, but in differ- 
ent manners such that combining them yields richer diagnos- 
tics then when using either line alone. This has not been done 
with the angular, spectral, and temporal resolution and cover- 
age furnished by the imaging spectroscopy presented here. In 
this initial paper we compare the chromospheric scenes seen in 
the two lines by displaying basic spectral measurements across 
a quiet-sun field of view. In particular, we show that both lines 
provide evidence of network heating, each in its own manner. 
So far, chromospheric heating in or above the magnetic con- 
centrations that cluster together to constitute the network has 
been most clearly evident as brightness enhancement of the cores 
of Can H&K. A recent discuss ion of its nature is given by 
iHasan & van BallegooiienI ( |2008|) . Additional information from 
other diagnostics is obviously desirable. 

There is an extended literature addressing chromospheric 
physics using filtergrams taken in Ha or Ca ii H & K, while most 
work using CaIR was performed with one-dimensional spec- 
trometry. We obtain a fresh view by combining Ha and CaIR 
in imaging spectroscopy, obtaining high angular resolution with 
adaptive optics at the Dunn Solar Telescope (DST) at the 
U.S. National Solar Observatory/Sacramento Peak and excellent 



spectral profile sampling by employing the Interferometric Bi- 
dimensional Spectrometer (IBIS) mounted at this telescope. 

High-resolution imaging spectr oscopy in Ca I R was ini- 
tiated with IBIS in the study of IVecchio etall (l2007l) and 
continued in the previous papers of this ser ies ( ICauzzi et al.l 
120081: IVecchio et al.ll2009l: iReardon et al.l l2009') mosdy address- 
ing the properties of the quiet chromosphere. Recent chro- 
mospheric studies with the Swedish 1-m Solar Telescope 
combine high cadence with 0.2 arcsec resolution but em- 
ployed Ha only in si ngle-wavelength filtergraph sampling (e.g 
let aLllr ""' - " ■ •■ 



Hansteen et al 112006; van Noort & Rouddc van der Voort 20()6|; 



De Pontieu et al.ll2007a ; .Rouppe van der Voort et al.ii2007i) . The 
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Dutch Open Telescope has registered man y Ha filtergram se- 
quences with subarcsecond resolution (e.g., iRutten et al.ll2008l) 
but only with coarse profile sampling and limited in duration 
by atmospheric seeing. Much longer multi-wavelength Ha fil- 
tergram sequences at this resolution were expected from the tun- 
able filtergraph onboard Hinode, but technical problems have 
severely limited its ability to sample the Ha profile at high 
cadence. The capability of IBIS and similar new Fabry-Perot 
instruments (the new Gottinge n FPI at the German Vacuum 
Tower Telescope described by iBello Gonza lez & Kneeij 20081 
and CRISP at the SST following the design of Schar me32006h to 
supply high-resolution imaging at fast cadence, with full profile 
sampling of both lines, and over fairly long durations represents 
a major step forward in investigations of the solar chromosphere. 

Our goal here is to show and compare what the quiet chro- 
mosphere looks like at such improved resolution and coverage 
of space, time, and wavelength in the two lines. Comparison be- 
tween the two is of interest because their formation differs in- 
trinsically, as discussed in Sect. 14. II Each has its advantages and 
disadvantages as a chromospheric diagnostic. We limit the anal- 
ysis here to a straightforward presentation of three basic profile 
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Fig. 1. Magnetic field evolution for the IBIS field of view (dashed circle), from GONG data at Ihttp : //gong . nso . edu" at 20- 
minute intervals. The time of observation is specified at the top of each panel. The IBIS run was from 17:50 to 18:40 UT. The 
contours show the apparent magnetic flux density at values of 10 and 20 Mx cm"^, for positive (solid) and negative (dotted) polarity. 



Table 1. IBIS profile scan parameters. 



Line 


Ha 


CaIR 


Nominal wavelength [nm] 


656.281 


854.214 


Instrumental profile FWHM [pm] 


2.2 


4.4 


Wavelength range [nm] 


0.21 


0.18 


Core sampling interval [nm] 


0.01 


0.008 


Number of wavelength steps 


22 


20 


Profile scan duration [s] 


5.2 


4.7 



measurements: the profile-minimum intensity, the Dopplershift 
of the profile mimimum, and the line-core width. Comparison of 
the general appearance of the two lines in these basic properties 
already yields informative similarities and differences. 



2. Observations, data reduction and profile 
measurements 

The imaging spectrometer IBIS is mounted at one of the two 
high-order adaptive-optics feeds at the DST. The latter, de- 
scribed by iRimme 13 (12004 . serves to stabilize the image and 
correct wavefront deformations in real time. IBIS consists of two 
Fabry-Perot interferometers in tandem and delivers sequences of 
narrowband images with the passband stepping quickl y through 
multip le se lectable spectral line s . IBIS is detailed in ICavallinil 
(120061) and lReardon & CavalUnil (|2008|) . 

In this paper we use data taken during 17:50-18:40 UT on 
15 March 2007. Table [T] specifies the spectral profile sampling. 
IBIS acquired 192 spectral image sets sampling the Ha, CaIR 
and Fe 1 709.0 nm lines with a cadence of 15.4 seconds between 
successive scans of all three lines. The Fei data are not con- 
sidered in this paper. Ha was scanned with a constant spectral 
sampling interval of 0.01 nm, CaIR with 0.008 nm sampling 
in the core and in each wing at three wavelengths separated by 
0.016 nm. 

The circular field of view has a diameter of 80 arcseconds 
and was located near disk center at heliocentric coordinates 
3.0 S, 10.1 W degrees, ju - 0.98, slightly eastward of a small 
bipolar plage. Figure [1] shows the magnetic content and evo- 
lution of the area from full-disk GONG magnetograms. It was 
very quiet, containing a few bipolar patches of magnetic con- 
centrations that appear as bright points in simultaneously taken 
G-band images. Only slight topology changes occurred during 
the acquired sequence. 



All images were corrected with appropriate dark and flat 
field calibrations. The transmission profile of the prefilter was 
measured and removed, resulting in good agreement between 
the spatio-temporal averages over the full sequence and the cor- 
responding profiles in the NSO/FTS solar spectrum atlas of 
[Wallace et al. ( 1998^ compare also the average CaIR profile in 
Cauzzi et al. 2008). 

In order to retrieve the line profile at each pixel, it is neces- 
sary to accurately align the images making up a spectral scan 
through a line, including "rubber-sheet" destretching to cor- 
rect remaining field-dependent deformations by seeing. The ref- 
erence for such alignment was generated from a sequence of 
broadband images obtained simultaneously with the narrowband 
images. The best broadband image every 15 seconds (i.e., during 
each full scan of all three lines) was found by measuring granular 
contrast. These 192 best broad-band images per spectral set were 
then co-aligned to each other, including destretching through lo- 
cal cross-correlation while taking care not to remove persistent 
solar flows. The destretched sequence was Fourier filtered in 
both space and time using a subsonic filter to remove the Fourier 
components corresponding to apparent horizontal motions faster 
than the photospheric sound speed (7 km s"'). This resulted in 
a high-quality master reference sequence which was then used 
to evaluate the seeing-induced local image motion, again using 
local-correlation tracking, for each individual broadband image. 
The destretch vectors that were derived for each of these were 
finally applied to the corresponding simultaneous narrowband 
image. This procedure results in stable, coaligned spectral-scan 
image sequences and so in reliable line profiles. During the mo- 
ments of best seeing the resulting image quality is close to the 
DST diffraction limit (0.22 arcsec for Ha, 0.28 arcsec for Ca IR). 

Each individual spectral profile at each pixel in each pro- 
file scan was interpolated onto a uniform wavelength grid with 
0.009 nm sampling for Ha and 0.008 nm sampling for Ca IR. 
The radially increasing instrumental blueshift, which results 
from the placement of the IBIS interferometers in a collimated 
mount, was also removed in this resampling. For each spec- 
tral profile we then determined the wavelength of the profile- 
minimum by fitting a second-order polynomial to the five spec- 
tral samples centered on the position with the smallest measured 
intensity. The minimum of this fit is taken as the instantaneous 
line center, yielding both its intensity (in relative data units) and 
its Dopplershift (in km s ' with respect to the profile-minimum 
wavelength of the spatio-temporal profile average over the flat- 
field data). 
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Fig. 2. Core width measurement. Left: Ha. 
Right: Ca IR. Eacli plot contains profiles from 
four different pixels to illustrate profile varia- 
tions. The thick grey curve is the solar atlas 
profile. The intensity scale is in relative units 
with respect to the atlas continuum value. For 
each profile /i(x, v, the wing intensities at 
given AAi from the profile minimum (marked 
by symbol pairs) are averaged. The core width 
is the wavelength separation of the profile sam- 
plings at half the intensity range between the 
profile minimum and this wing average. The 
same pixels are selected in Fig. [7] 



The measurement of line width is less straightforward. The 
scan extents specified in Table [T]do not reach the local continua 
outside the lines, so that we cannot measure the profile width 
at a fixed fractional depth from the continuum. Furthermore, the 
shape of these chromospheric lines often varies significantly, es- 
pecially in their cores. For Ca IR in particular, the extended pho- 
tospheric wings tend to remain relatively stable while the chro- 
mospheric core varies dramatically. Expressing the latter's chro- 
mospheric properties in terms of photospheric continuum inten- 
sity thus is not a desirable approach. We therefore use only the 
inner, chromospheric part of the profile to define the property we 
call core width. 

We evaluated these core widths per instantaneous profile per 
pixel by first determining and averaging the wing intensities at 
a specified wavelength separation A/l/ from the measured wave- 
length of the profile minimum. We then defined and measured 
the core width as the wavelength separation of the two profile 
flanks at half the intensity range between the minimum and this 
wing intensity. The separation parameter was set at AAi - 90 pm 
for Ha and AAi = 60 pm for CalR. These values were selected 
through inspecting many profiles with the aim to separate the 
chromospheric core from the photospheric wings as well as 
possible. The choice remains somewhat arbitrary, but we feel 
that this procedure is a better one than halfwidth measurements 
higher up in the line that would mix disparate photospheric and 
chromospheric contributions. Figure [2] iflustrates the procedure 
for the Ha and CaIR profiles from four selected pixels at one 
moment. Note that the CaIR core varies much more than the 
Ho- core, and that the highest Ca IR profile has a rather unusual 
redshift and blue-wing brightening to which we return below. 

Finally, the observed intensities were converted into bright- 
ness temperatures by equating the spatio-temporal average of 
each line minimum to the corre sponding absol ute intensity value 
in the solar spectrum atlas of Neckel (1999) which is a cali- 
brated version of the atlas of Wallace et al. (1998i) , and then ap- 
plying the inverse of the Planck function. These formal tempera- 
tures represent kinetic temperatures only when LTE holds for the 
source function, which is not the case for either line (Fig.[8]l. The 
conversion serves primarily to undo the Planck-function sensi- 
tivity variation with wavelength. 



3. Parameter displays and comparisons 

Figure |3] shows our measurements from a single profile scan of 
the two lines, taken at 17:52:34 UT during one of the best-seeing 
moments. For each line the measurements are as specified in the 



previous section: the intensity of the line profile minimum ex- 
pressed as brightness temperature, its Dopplershift, and the line 
core width defined as in Fig. |2] We describe them column by 
column, comparing each panel with pertinent other ones. These 
comparisons are quantified in the form of pixel-by-pixel corre- 
lation diagrams in the lefthand part of Fig. |6] which we describe 
in parallel. The diagrams are constructed from 66 spectral sets 
taken during very good seeing to obtain the best statistics. The 
caption explains the format, while the axes provide the numer- 
ical ranges of the greyscale variations in Fig. |3] The innermost 
contours of the diagrams in Fig. |6] represent the maxima of the 
overall distribution, while the curves on the top and right of each 
figure show the distribution of the relative quantity on each axis. 

The two images of the profile-minimum intensity at left in 
Fig. [3] differ dramatically in their portrayal of the solar atmo- 
sphere. The corresponding scatter diagram (first panel of the 
second row in Fig. |6]l indeed shows a total lack of brightness 
correlation between the two, except for the very brightest pix- 
els. The Ha image shows masses of bright and dark chromo- 
spheric fibrils, sometimes with marked brightening towards one 
end. Comparison with the magnetograms in Fig. [T| shows that 
these bright endings correspond to the stronger magnetic con- 
centrations, but the magnetic network pattern cannot be recog- 
nized directly from the Ha image. That is much more easily 
done in the CaIR intensity image, in which the extended dif- 
fuse bright patches constituting the chromospheric network cor- 
respond reasonably well to the unsigned magnetogram distribu- 
tion. The most active patch is located at the field edge at about 
2 o'clock; it contributes most of the slight bright-bright correla- 
tion in Fig. |6] with brightness temperature T\, x 4300 K in both 
lines. The dark internetwork areas contain small-scale brighten- 
ings that la st only a min ute or so and mark acousti c shocks (see 
Movie 1 in lCauzzi et"aI.L2 008: Vecchio et al.ll2009h . 

We now turn to the Dopplershift maps in the center pan- 
els. The Ha Dopplermap has large similarity to the Ha profile- 
minimum map in its overall spatial structure. It shows the 
same masses of fibrils, roughly with the same bundle patterns 
and orientations. However, the two are dissimilar in detail: the 
Dopplermap fibrils tend to be shorter and they are dark or bright 
(moving upward or downward) without obvious relation to the 
contrast in the profile-minimum image, except for the bright- 
est and the darkest features whose co-spatiality indicates ex- 
ceptional downdrafts in fibril ends near the network. The lack 
of pixel-by-pixel similarity is quantified in the second panel of 
Fig.|6]which indeed shows only very slight correlation. The cor- 
responding CaIR diagram in the first panel of Fig. |6] shows 
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Fig. 3. Image representations of profile measurements at a moment of good seeing (scan nr 11). Solar North is to the top. West to 
the right. The greyscale is bytescaled separately for each panel, for the Dopplershifts symmetrically with respect to zero (byte value 
127; blueshift black, redshift white). Upper row. Hcc. Lower row: CalR. First column: minimum intensity per pixel, measured as 
brightness temperature. Second column: Dopplershift of the profile minimum per pixel. Third column: core width, measured as in 
Fig. El 



slightly higher coherence but with negative correlation: the in- 
tensity tends to be higher at locations with larger blueshift (dark 
in the Dopplermap), except for the brightest pixels which there- 
fore diminish the negative correlation in the overall Pearson co- 
efficient. We show below that the countercorrelation comes from 
the internetwork shocks. 

Much better coherence is shown in the second panel of the 
second row of Fig. |6] which plots the per-pixel correlation be- 
tween the two Dopplershift measures. One would not expect 
such high correlation comparing the center-column images in 
Fig. |3] since they appear rather dissimilar, but closer inspec- 
tion shows that indeed the CaIR Dopplermap has similar fib- 
rils around network as the Ha Dopplermap, although less ex- 
tended. In the internetwork the former shows many roundish 
patches that must correlate fairly well with the more fibrilar- 
shaped Dopplershift features in Ho- to produce the good over- 
all correlation in the scatter plot. This correlation may be pro- 
duced by acoustic shocks seen simultaneously in bo th lines, or if 
shocks from below buffet Wa fibrils as suggested by Rutten et al.j 
( 2008 ) or excite local turbulence as suggested by Reardon et al,: 
(l2008h . The two Dopplermaps so possess larger pixel-by-pixel 
correlation than cursory inspection suggests, just the opposite of 
the Ho- brightness-Dopplermap comparison. 

The remaining two panels in Fig. [3] and the remaining four 
panels in the lefthand part of Fig. |6] concern the core width mea- 



sures of the two lines. It is most striking that the Wa width pat- 
tern has much larger similarity to the Ca IR minimum intensity 
than to anything else; this is the key point of this paper. They 
both show the network as bright patches on a dark background, 
with longer fibrilar extensions in Ho- width than in Ca IR bright- 
ness. The internetwork grains stand out more clearly in CaIR, 
but there is very good similarity in the underlying internetwork 
patterns. This large apparent correspondence is confirmed on the 
pixel-by-pixel level by the second scatter diagram in the bottom 
row of Fig.|6]which shows both bright-bright and dark-dark cor- 
relation with high significance. There is also significant correla- 
tion between Ca IR brightness and width and between Ha width 
and Ca IR width in the first column of Fig. |6] There is no corre- 
lation at all between Ha intensity and width, in accordance with 
the absence of correlation between the two brightness images. 

Some peculiar patches appear at times in the Ca IR parameter 
images. A clear example is given by the feature near (32", 73") 
in Fig. |3] which appears very bright in the minimum intensity 
and in the line-width panels, while being very dark (upflow) with 
a surprising bright center (downflow) in the Dopplershift map. 
These features often represent particularly strong internetwork 
shocks with a strongly asymmetric CaIR line core profile for 
which our simple measurements are inadequate. However, such 
occurrences are rare (much less than one percent), and do not 
affect the statistical correlations in Fig.|6] 
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Fig. 4. Temporal averages over the 48-min sequence duration for the central part of the field of view, excluding moments of 
bad seeing (about one-third of the profile scans). The layout furnishes the same quantities as in Fig. [3] Each panel is bytescaled 
separately, the Dopplershifts symmetrically with respect to zero. The two markers along the y-axes specify the cut locations used in 
Fig. |5] They pass through a patch of network area in the center and a weaker field concentration near x - 63 arcsec, and through 
two dark patches of small Ca IR core width in the last panel. The four markers along the x- axes specify pairs of pixels along each 
cut for which the spectral profile evolution is shown in Fig.|2l 



Figure |4] shows temporal averages over the sequence dura- 
tion of the same quantities as in Fig. [3] Moments of less good 
seeing were excluded by setting a quality threshold at the mean 
rms variation of the Ho' profile-minimum intensity images, dis- 
carding nearly one-third of the samples, spread fairly evenly 
in time. The righthand half of Fig. |6] contains the correspond- 
ing pixel-by-pixel scatter plots for the whole IBIS field of view, 
while only the central part of the field is shown in Fig. ID This 
is done to enhance the spatial scale in Fig. |5j which shows the 
temporal evolution for the pixels along the two fixed-y cuts in- 
dicated by the horizontal white markers in Fig. ID They are se- 
lected to pass through the two distinct black patches in the last 
panel of Fig.|4]and are 1 .5 arcsec apart. We describe these figures 
together. 

Figure |5] shows that periodic phenomena abound. All time 
slices show ubiquitous oscillatory modulation, not only in the 
internetwork areas but also for the fibrils that jut out from the 
small area of active network sampled in the middle. The oscilla- 
tory signals are most obvious in the Dopplershift patterns (center 
column) which are remarkably similar for the two lines. They are 
least clear but still visible in Ha intensity (first panel). They aver- 
age out when taking a temporal mean, so that the spatial patterns 
defined by longer-lived structures and motions gain visibility in 
Fig.|4]in comparison with Fig. [3] 



In all measures the internetwork is full of three-minute os- 
cillations with shock signatures, cl earest in the CaIR i ntensity 
panels in the form of bright grains jVecchio et alj|2009l) . These 
are comparab le to the Can H2V and K21/ gr ains reviewed by 
Rutten & Ui tenbroekl d 1 99 lb and explained bv lCarlsson & SteirJ 
( I l997h as steepening and vertically interacting shock waves that 
travel up until wave conversion occurs at the magnetic canopy 
dBogdan et al. I l2003h . The Ha time slices show that Ha par- 
takes in this three-minute modulation, even in its intensity al- 
though that has no resemblance nor correlation to CaIR in 
the instantaneous images and scatter diagrams. Further study 
of this modulation calls for Fourier analysis including phase- 
differe nce evaluation with n etwork/internetwork separation fol- 
lowing iKrijger et al.l ( 1200 ll) . which we postpone to a future pa- 
per 

The network patch at the center of the time slices shows os- 
cillatory behavior of its own, most clearly in the Dopplershift 
panels. Some fibrils that jut out from it, especially towards the 
left in Ha, have alte rnating flows and may repre sent dynamic fib- 
rils as d escribed by iHansteen et al.l (l2006l) and iDe Pontieu et al.l 
(l2007al) . The central part d isplays slower undulations as de- 
scribed bv lLites et al.l (Il993l) . The fibril-dominated regions sur- 
rounding the network show a distinctly different behavior than 
either the nework or internetwork regions, as described by 
IVecchio etai] (120091) . 
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Fig. 5. Space-time slices showing temporal behavior along the two cuts in the x-direction that are defined in Fig. |4]by white 
horizontal markers along the y-axes. The upper half of this figure corresponds to the higher cut. The six panels per cut show the 
same quantities as in Figs.[3]and|4l Each panel is bytescaled separately, the Dopplershifts symmetrically with respect to zero. The 
white markers along the spatial axes specify four pixels, two per slice, for which the spectral profile evolution is shown in Fig.|2] 
The white markers near the bottom of the time axes specify the moment of observation of the profiles in Fig. |2] and the parameter 
displays in Fig. [3] 



The temporal averaging in Fig. |4] and the righthand part of 
Fig. |6] reduces the spread due to these oscillations and other 
short-lived phenomena, and so serves to enhance longer-term 
structuring and relationships (note also the reduced ranges of the 
axes in the right-hand side of Fig.|6]l. The averaging over the os- 
cillatory components indeed reduces the blue-bright correlation 
between Ca IR brightness and Dopplershift in the first panel of 
Fig.|6l introduced by internetwork acoustic shocks. In the tempo- 
ral averaging of the Ho- data, the fibril contribution gains domi- 
nance and improves the correlation with intensity, corresponding 
to larger resemblance between the first two Ha panels in Fig.|4] 
than in Fig. [3] 

The scatter plots in the second and third rows of Fig. |6]do 
not change much in the temporal averaging. The Dopplershifts 
comparison in the second row looses the internetwork oscilla- 
tions but maintains fairly good correlation, indicating a corre- 
spondence in the fibril flows also seen in the center panels of 
Fig- El In the third row there is slight improvement in the high- 
high correlation between Ca IR brightness and width. 



The most striking change that the temporal averaging brings 
is increased similarity between the panels relative to the Ha 
width, CaIR profile-minimum intensity and CaIR width of 
Fig.|4]in comparison with Fig. [3] In particular, the time-averaged 
Ha width (third panel in upper row) shows a pattern of smooth, 
bright network aureoles on a dark internetwork background that 
appears very similar to the Ca IR intensity (first panel in lower 
row). The bright network patches seem smaller in the latter, but 
this results from the non-linearity displayed in the last panel of 
Fig.|6] which quantifies this comparison. It shows high correla- 
tion along the whole scatter-sample mountain, but the turnover 
of the upward tail implies that the brightest network patches ap- 
pear larger in Ha width. The whole mountain ridge is apprecia- 
bly narrower in the time-averaged correlation (last panel) than in 
the instantaneous one (bottom second column). Thus, both net- 
work (contributing the upper part) and internetwork (lower part) 
gain pattern similarity from temporal averaging in this compari- 
son. This tight relationship suggests a common physical origin, 
which we identify as temperature below. The width-width cor- 
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Fig. 7. Spectrum- versus-time plots showing the spectral behav- 
ior of Ha (upper row) and Ca IR (lower row) for the four se- 
lected pixels specified by markers in Figs. |4]and|5] The wave- 
lengths along the abscissae are expressed as Dopplershift from 
the nominal line center. The logarithm of the spectral brightness 
is shown with the same greyscale for all panels along each row; 
it is logarithmic in order to detail line-core variations, especially 
for Ha which is too uniformly dark otherwise. The first and sec- 
ond columns are for pixels along the lower cut in Figs.|4]and|5] 
The first pixel samples normal network, the second the network 
feature with low Ca IR core width. The third and fourth columns 
are for pixels along the upper cut and sample the near-network 
feature with low Ca IR core width and internetwork further away. 
The black markers near the bottom of each panel specify the time 
of observation selected in Figs. |2] and [3] The second pixel hap- 
pened to lie in a dynamic fibril at that moment, causing large 
Ca IR redshift and blue-wing brightening also seen in Fig.|2] 



relation (bottom panel of the third column in Fig. |6]l is just as 
good as the Ha-width vs. Call intensity correlation for large 
core widths but it vanishes for small core widths, indicating that 
temperature dominates the Ca IR width only when high and that 
temporal averaging over the internetwork shocks weakens their 
temperature signature, perhaps through their short duration with 
respect to their periodicity. 

Note that the Ca IR width image in Fig.|4]shows dark centers 
in some of the bright network patches. The yet brighter network 
patches near 2 o'clock and 6 o'clock at the edge of the IBIS 
field of view, not included in Fig.|4] have similar dark cores. The 
lower slice in Fig. |5] samples the darkest network core in Fig.|4] 
The upper slice cuts through a comparably dark feature in Ca IR 
width that lies adjacent to the network patch. However, com- 
parison of the two dark Ca IR width patches in the time slices 



in Fig. |5] (the two last panels of each half) suggests that they 
have intrinsically different nature. The network core patch in the 
lower panel lies at the center of the fibril rosette and is partially 
visible also in Ha width. The dark patch adjacent to the network 
in the upper panel instead seems to be regularly oscillating inter- 
network, only darker overall. 

Figure |7] adds further information on these features by com- 
paring the spectral behavior of the two lines with time, respec- 
tively for pixels within these two dark patches (center columns) 
and for comparable pixels elsewhere. The two left-most columns 
sample network pixels, the two right-most ones sample internet- 
work. The first thing to note about the Ha line core is its much 
larger width with respect to Ca IR, and its ubiquitous presence 
(note also that the spectral scale of Fig.|2]is different for the two 
lines, covering about 0.2 nm around Ha and 0.18 nm around 
CaIR). This basal core width is widest in the first panel, while 
the superimposed Doppler excursions become a major compo- 
nent only in the internetwork (rightmost columns). The CaIR 
core shows these Doppler excursions as well, but more visibly 
due to its reduced width. Regular shock trains occur most clearly 
in the internetwork panels for the Ca IR line and are also rec- 
ognizable in the corresponding Ha panels. In both dark CalR- 
width patches (second and third columns) the Ca IR core indeed 
appears generally narrower, but also less dark than for the com- 
parison pixels in the first and last columns. Profile inspection 
along the second column shows that in these pixels CaIR of- 
ten has a very shallow chromospheric core for which our width 
algorithm produces small values. 

Finally, we note that the time of observation selected for 
Fig.|2]happened to sample a feature indicative of dynamic fibrils 
in the second panel of Fig. |7] i.e., displaying large blue-to-red 
successive shifting in Ca IR. It has correspondingly large redshift 
and blue-wing brightening in Fig.|2] There is no clear signature 
of this feature in Ha. 



4. Discussion 

4.1. Temperature sensitivity 

The Ca IR and Ha lines differ in two significant ways in their 
temperature sensitivity. First, Ca IR is much closer to LTE both 
with regards to opacity and to source function behavior. The lat- 
ter difference is demonstrated in Fig.[8]which represents a didac- 
tic step beyond the LTE comparison in Fig. 6 of L eenaarts et al.l 
(120061) . by comparing standard NLTE formation for the two 
lines. It shows dramatic difference between the two lines in 
their sensitivity to the temperature in the outer atmosphere. In 
ea ch case the lines are computed in NLTE with the RH code 
of lUitenbroe3(l200lh . using a sufficiently complete model atom 
and accounting for coherent scattering in Ly a, Ly y6 and H&K 
while assuming complete redistribution for all other lines. In the 
upper panel, for the KURUCZ radiative-equilibrium model Ha 
is a deep-photosphere line of which the source function decou- 
ples from the Planck function as deep as log rsoo ~ 0, together 
with the Balmer continuum. The latter. Ha, and Lya combine 
together in setting the pertinent hydrogen populations through 
intricate interactions involving Balmer continuum pumping. Ha 
photon l osses, Rydberg ladder flo w, and detailed balancing in 
Lya (see lRutten & Carlssonll 1994b . The Ha source function S ,i 
closely follows the mean intensity J,i (not shown) as if it were 
an ordinary two-level scattering line. Hence, its temperature sen- 
sitivity responds mostly to deep-photosphere perturbations. In 
contrast, the Ca IR line has much larger opacity, so that its source 
function follows the temperature out to log rgoo ~ -4 before the 
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Fig. 8. Demonstration of difference in temperature sensitivity 
between CaIR and Ha. Each panel shows the temperature strat- 
ification of a standard solar model atmosphere and the resulting 
total source functions 5^ at the nominal line-center wavelengths 
for CaIR and Ha, as function of the continuum optical depth 
at /I = 500 nm and with the source functions expressed as for- 
mal temperatures through Planck function inversion. The arrows 
mark r = 1 locations. Upper panel: radiative-equilibrium model 
KURUCZ from Kurucz (1979, 1992a, 1992b). It was extended 
outward assuming constant temperature in order to reach the op- 
tically thin regime in Ca IR. Lower panel: em pirical continuum- 
fitting model FALC of lFontenla et al.l ( Il993h . Its very steep tran- 
sition region lies beyond the top of the panel but causes the near- 
vertical source function increases at left. 
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Fig. 9. Width- width correlation between Ha and CaIR. The 
contours reproduce every second one in the next-to-last panel 
of Fig. |6] The dashed curve represents Eq. (HJ with best-fit pa- 
rameters. 



onset of its scattering drop with 5^ ~ J,\ which makes the line 
very dark. 

In the lower panel the photospheric part of the FALC temper- 
ature stratification is closely the same as the KURUCZ one and 
so the formation of the lines there remains similar, with thermal- 
ization out to log tsoo ~ for Ha and log rgoo ~ -4 for Ca IR. The 
addition of the FALC chromosphere results in a source function 
bulge for CaIR that is familiar from the older explanations of 
Ca II H (fe K reversals before partial redi stribution and shock dy - 
namics were included (e.g.. Fig. 13 of iLinskv & AvretBll970l) . 
For Ha, the chromosphere plus transition region act as the addi- 
tion of a barely opaque, hot, radially inhomogeneous slab well 
above the photosphere. It raises the scattering source function 
across the upper photosphere but without giving it sensitivity to 
the local temperature since there is virtually no local Ha opacity 
there. 

The T = 1 arrows in Fig. |8] indicate that the predicted line- 
center brightness temperatures ~ 7'rad(T"^ = 1) happen to 
be nearly equal for Ha at ~ 3900 K for both models, but 
in very different manner. The photosphere-only prediction for 
Ca IR line-center brightness (Tb ~ 2500 K) is instead much lower 
than with a chromosphere (Tb ~ 3600 K). The observed mean 
values from Fig.|6]are Tb -4100 K for Ha and ~ 3800 K for 
Ca IR. Given the very different structure of the source function 
of the two lines, the closeness of their two brightness tempera- 
tures in the model with a chromosphere is more likely a coin- 
cidence imposed by the atmospheric model than proof that the 
lines supply temperature information from the same height in the 
atmosphere. The Ha line-minimum intensity is much more sen- 
sitive to deep-photosphere temperature perturbations, whereas 
the Ca IR line minimum has more sensitivity to chromospheric 
temperature perturbations. The first columns of Figs. [3] and |4] 
confirm this difference dramatically by displaying very dissimi- 
lar chromospheric intensity scenes, not only in the instantaneous 
comparison (Fig. [3j but also in the time-averaged one (Fig. |4]i. 
In contrast, the closeness of the formation heights in the chro- 
mospheric model suggests that the Dopplershifts and Doppler 
broadening, which are encoded locally at the last photon scatter- 
ing independent of the source function nature, should be much 
more similar. This may explain the larger similarity between the 
two lines in the corresponding panels of Figs.|3]-|6] However, the 
large difference in Ha formation height between the KURUCZ 
and FALC models on the one hand and the larger apparent extent 
of fibrils in Ha than in Ca IR in Figs. [3]and|4]on the other hand 
imply that there is no guarantee of Doppler-encoding similarity 
between the two lines. 

The second major difference between the two lines is that 
they differ by a factor of 40 in atomic mass. This affects the 
thermal Doppler broadening and therefore the balance between 
thermal and nonthermal Doppler broadening of the line profile. 
The latter was traditionally described by the ad-hoc fudge pa- 
rameters microturbulence and macroturbulence, with the former 
entering the Dopplerwidth and the latter applied as Gaussian 
smearing of the computed emergent line profile. In recent times 
the usage of these parameters has been replaced by structure- 
resolving line synthesis from numerical simulations, but so far 
only for granular convection and magnetoconvection in the low 
photosphere. Given the low atomic mass. Ha is likely to have 
thermal broadening dominate over systematic motions, whereas 
the reverse holds for Ca IR. The spectrum-versus-time displays 
in Fig. Qconfirm this difference qualitatively. Hence, we attribute 
the large, generally-present core width of Ha to thermal broad- 
ening, and therefore the bright patches in the third panel of Fig.|4] 
to high temperatures which set the core width locally, at the last 
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scattering. The good correspondence of the Hor hne width with 
the patches of high Ca IR profile-minimum intensity in the fourth 
panel then results from the appreciable CaIR source function 
sensitivity to temperature. 

Figure |9] quantifies this conclusion in the form of a simple 
model which reproduces the observed correlation in the next-to- 
last panel of Fig. |6] The dotted line in Figure plots the widths of 
the two lines against each other that were computed from 



AA 



w 



+ . + vr . 

micro intnns 



(1) 



for temperature T increasing from r^in = 5 000 K to r^ax - 
60 000 K, a thermal contribution Vtherm - yJlkT/m where k is 
the Boltzmann constant and m the atomic mass, a microturbu- 
lent contribution Vy^icm = 1 + 10 {T -T^,n)l(T^^^-T^^„) km s"', 
and an intrinsic contribution Vinmns set to 43 km s"' for Yia and 
16 km s ' for CalR. The latter two contributions were adapted 
to fit the observed correlation. The temperature-dependent mi- 
croturbulence increases from 1 km s ' to 11 km s in rough 
agree ment with older estimations from ultraviolet line broad- 
ening ('Canfield & Beckers' ' 19761: iTripp et alJI 19781: cf. Fig. 11 
of [Vernazza et al. 1981). The intrinsic widths represent base 
contributions making up for the neglect of radiative transfer 
Eq. ([T]i may be seen as describing line broadening in a Schuster- 
Schwarzschild slab or cloud with the irradiative profile and in- 
ternal line transfer mimicked by these intrinsic widths. It is only 
a simple ad-hoc fit, but it demonstrates that the observed mean 
width-width correlation in Fig.|6]may be attributable to temper- 
ature. The value Tmax - 60 000 K seems very high in comparison 
with standard one-dimensional models in which both hydrogen 
and once-ionized calcium are fu lly ionized at s uch temperatures, 
but observations as in Fig. 13of lRuttenl(l2007h suggest the pres- 
ence of Ha even at very high temperature. 

Note that such core width behavior cannot be reproduced 
by specification of a deeper temperature rise as done by 
iFontenla et al] ( Il993l) . Their A, C, . . . , F and P models were de- 
signed to represent progressively hotter chromospheres in quiet 
cell interior, average quiet sun, network, and plage by placing 
the chromospheric temperature rise at increasingly larger col- 
umn mass. However, when measuring with our method the line 
widths from profiles synthesized in NLTE (as in Fig. [8]), this se- 
quence produces smaller core width with increasing activity in- 
stead of the actual core-width increase we observe. 



4.2. Chromospheric heating 

The demonstrations above make us conclude that the high cor- 
relation between time-averaged \ia core width and CaIR line- 
minimum intensity represents common sensitivity to tempera- 
ture. Hence, the network patches display excess heat. The gen- 
eral, basal presence of a wide Ha core in the first columns of 
Fig. [T] as well as the large Ha widths ubiquitously present in 
the network in Fig. |5] indicate that this heating is rather constant 
in time in the magnetic concentrations. Temporal averaging over 
episodes of strong heating occurring at the base of the fibrils jut- 
ting out from the network (cf. the rightmost panels of Fig. |5]l 
produces the smooth bright network aureoles in the third panel 
of Fig.|4] The dark network patch of small CaIR width sampled 
in Figs. |5] and Fig. |7] indicate reduce core-widths related to a 
significant weakening of the line core, likely marking maximum 
heating at the center of the network patch. 

The calculated fit to the width-width correlation in Fig. |9] 
includes a sizable microturbulence component which scales 



with temperature. Its necessity suggests the action of unre- 
solved waves. Perhaps these are the ones affecting or produc- 
ing the "straws" or "sp icules-I I " obs erved in Can H, respec- 
tively near the limb bvlRuttenI (l2006h and across the limb by 
iDe Pontieu et al. ( 2007bl). th at were interpreted as Alfven waves 



bv lDe Pontieu et aTl (l2007c). They may also exp lain unresolved 
ultraviolet line broadening ( Mcintosh et al.ll2008h and the "rapid 
bl ueshift events" at t he edg e of network concentrations reported 
bv lLangangen et all (l2008l) . It is also interesting to note that the 
higher-temperature regions (i.e. network and fibrils) display a 
stronger tur bulence spectrum of the Ca IR velocity than the in- 
ternetwork jReardon et al.ll2008l) . 

The similarity between Ha width and CaIR intensity ex- 
tends to the oscillation sequences in the internetwork areas in 
Fig. |5] The Ha width shows internetwork grains less markedly 
than CaIR intensity, but it follows the same modulation. The 
bright phases may therefore be interpreted as localized temper- 
ature increases, but the overall darkness of the internetwork in 
Fig. ID implies less total heating than in the network. The dark, 
near-network patch in CaIR width sampled in Figs.|5]and Fig. [7] 
also seems to have smaller CaIR width from being brighter 
than the comparison internetwork pixel, perhaps from more sub- 
canopy heating so close to network. 

4.3. Future interpretation 

The next step in profile interpretation is to distill more in- 
trinsic line formation parameters. Traditionally, chromospheric 
fine structu re is d escribed through cloud modeling (see re- 
view by ' Tziotzioul i 2007). In its most basic form it delivers the 
four parameters: optical thickness, source function (sometimes 
parabolic), Dopplershift, and Dopplerwidth for a homogeneous 
cloud irradiated from below. Our finding that in particular Ha 
profile-minimum intensity has no correlation with the actual 
temperature (as sampled by Ha width) nor with the instanta- 
neous Ha Dopplershift implies that the rich fibrilar structure 
displayed by this line in narrow-band filtergrams is indeed a 
complex mixture of all four parameters. A complication adding 
uncertainty is the required specification of the Ha profile in 
the cloud irradiation from below and aside. Figure [8] suggests 
that combination of radiative-equilibrium modeling with back- 
radiation from the cloud may serve as first approximation. 

Granted cloud-model determinations of these four parame- 
ters, the next issue is how to interpret the cloud source func- 
tion. Already in his original formulation for off-limb spicules, 
Beckers (e.g., 1964, 1968, 1972) converted the basic empirical 
cloud parameters i nto de nsity and temperature through the ta- 
bles of iGiovanellil ( Il967h for NLTE hydrog en ionization. More 
recen t reformulations of this approach (e.g.. lAlissandrakis et al.l 
1 19901 : iTsiropoula & Schmiedeiill997l) have generally been re- 
placed by large-volume NLTE line profile computation followed 
by data "invers ion" through best-fit par ameter adaptation (see 
lTziotzioul2007l for references). However. iLeenaarts et al.l(l2007l) 
have undermined this technique by finding that in chromospheric 
structures which suffer repetitive episodic heating, the actual 
populations of the atomic levels governing Ha can be very 
far from time-independent statistical equilibrium. Since shocks 
seem ubiquitously present in the chromosphere (center column 
of Fig. |5]l. Ha cloud-model interpretation needs to be reformu- 
lated including time-dependent hydrogen balancing. A first ap- 
proximation may be to maintain the high degree of ionization 
and the corresponding large Ha population fractions reached at 
moments of high temperature. 
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Obviously, it will be even better to compare data as ours di- 
rectly to numerical spectral line synthesis employing ab-initio 
MHD simulations of time-dependent chromospheric fine struc- 
ture, but these are not yet available (apart from unrealistic ones 
assuming LTE even in the chromosphere). When realistic ones 
with proper Ca IR and Ha synthesis arise, reproducing the scat- 
ter diagrams in Fig.|6]will provide an excellent proving ground. 

5. Conclusion 

Imaging spectroscopy with IBIS provides diagnostics well be- 
yond slit-limited spectroscopy or bandpass-limited filter imag- 
ing. The long-duration simultaneous CaIR and Ha profile- 
resolved IBIS sequences presented here give direct evidence of 
chromospheric heating in network. The clincher measurement is 
the Ho- core width. 

The network heating appears to occur fairly continuously 
and to affect also the nearby base of fibrilar features jutting out 
from it, but it does not extend as far as, nor shows the fibrilar fine 
structure of, the petals of classical rosettes in Ha filtergrams. 
Their visibility arises from complex combinations of Doppler 
width, Dopplershift, optical thickness, and source function vari- 
ations, without direct coupling to fibril temperature. In contrast, 
the striking correspondence between Ha core width and Ca IR 
line-minimum brightness (Figs.|4]and|6]l proves that the latter is 
a better temperature proxy. The best, however, is the Ha core 
width. 
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Fig. 6. Scatter-plots between measured quantities for all pixels in the full field of view. The first two columns are instantaneous 
pixel-by-pixel comparisons between the quantities specified along the axes, for the 66 profile scans with the best seeing (10^ 
samples). The third and fourth columns are similar comparisons for the same quantities but temporally averaged over the whole 
sequence, excluding the third of the scans with the worst seeing (1.6 x 10^ samples). The first two columns correspond to Fig. [3] the 
other two to Fig.|4]but sampling the full field. Plot format: to avoid plot saturation, single pixel-by-pixel sample plotting is replaced 
by plotting sample density contours, except that the extreme individual pixel-by-pixel comparisons beyond the outer contour are 
shown in columns 3 and 4. Extended arcs in these outer scatter clouds come from isolated, non-characteristic features. The solid 
curves along the top and right sides show the occurrence distribution per quantity. The dashed curves specify the first moments 
along cuts in each axis direction. They align at perfect correlation (tilt to the upper right) or anticorrelation (tilt to the upper left) and 
they are perpendicular along the axis directions in the absence of any correlation. The numbers near the lower-right corners specify 
the overall Pearson correlation coefficient. The intensities of the profile minima are quantified as brightness temperatures. Positive 
Dopplershift implies downdraft. The scales differ between the instantaneous and sequence-averaged comparisons, especially for 
Dopplershift. The contours in the latter are noisier due to the smaller statistics. Top row: the intensity of the profile minimum 
against its Dopplershift for CaIR and Ha. Second row: intensity-intensity and Dopplershift-Dopplershift comparisons. Third row: 
minimum intensity against core width for CaIR and Ha. Bottom row: Ho- core width against CaIR core width and against CaIR 
minimum intensity. 



